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SUMMARY 


This paper treats the kernel function of an integral equation that 
relates a known or prescribed downwash distribution to an unknown lift 
distribution for a harmonically oscillating finite wing in compressible 
subsonic flow. The kernel function is reduced to a form that can be 
accurately evaluated by separating the kernel function into two parts: 
a part in which the singularities are isolated and analytically expressed 
and a nonsingular part which may be tabulated. The form of the kernel 
function for the sonic case (Mach number of l) is treated separately. 

In addition, results for the special cases of Mach number of 0 (incom- 
pressible case) and frequency of 0 (steady case )' are given. 

The derivation of the integral equation which involves this kernel 
function, originally performed elsewhere (see, for example, NACA Technical 
Memorandum 979) is reproduced as an appendix. A second appendix gives 
the reduction of the form of the kernel function obtained herein for the 
three-dimensional case to a known result of Possio for two-dimensional 
flow. 


INTRODUCTION 


The analytical determination of air forces on oscillating wings in 
subsonic flow has been a continuing problem for the past 30 years. 
Throughout the first and greater part of this time, efforts were directed, 
mainly toward the determination of forces on wings in incompressible flow. 
These efforts have led to important closed-form solutions for rigid wings 
in two-dimensional flow (ref. l), to solutions in terms of series of 
Legendre functions for distorting wings of circular plan form (refs. 2 
and 5)> and to many approximate, yet useful, results for wings of elliptic, 
rectangular, and triangular plan form (see, for example, refs. 4 to 12). 
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Although these results for incompressible flow play a highly signif- 
icant roll in applications of unsteady aerodynamic theory, the advent of 
higher and higher speed aircraft during the last 15 years has brought a 
growing need for knowledge of the effect that the compressibility of air 
might have on unsteady air forces, or for analytically derived unsteady 
air forces based on a compressible medium. The transition to results 
for a compressible fluid from those for an incompressible fluid is not 
likely to be accomplished by applications of simple transformations or 
correction factors, such as the well-known "Prandtl-Glauert factor for 
steady flow. This difficulty is associated with the fact that the time 
required for signals arising at one point in the medium to reach other 
points gives rise not only to changes in magnitudes of forces but also to 
additional phase lags between instantaneous positions, velocities, and 
accelerations of the wing and the corresponding instantaneous forces 
associated with these quantities. In order to obtain results for the 
compressible case, it therefore appears necessary to deal directly with 
the boundary -value problem for this case. 

The boundary -value problem for a two-dimensional wing in compressible 
flow has been successfully attacked from two points of view. First, by 
consideration of an acceleration or pressure potential, Possio (ref. 13 ) 
reduced the problem to that of an integral equation relating a prescribed 
downwash distribution to an unknown lift distribution. The kernel of 
this integral equation, which is a rather abstruse function, was reduced 
to a form that, except at singular points, could be evaluated. Schwarz 
(ref. Ik) later isolated and determined the analytic behavior of the 
singular points of Possio^ results and made fairly extensive tables of 
the kernel function. These tabular values were used by various investi- 
gators (for examples, refs. 15 and l6) to obtain, by numerical procedures, 
initial tables of force and moment coefficients for oscillating wings in 
compressible subsonic flow. 

The second successful approach to the solution of the boundary -value 
problem for a two-dimensional wing (see refs. 17 to 19) is achieved by a 
transformation to elliptic coordinates followed by a separation of vari- 
ables that reduces the boundary -value problem from one in partial- 
differential equations to one in ordinary differential equations of the 
Mathieu type. The solutions turn out as infinite series in terms of 
Mathieu functions. Numerical results obtained recently by this proce- 
dure agree with results previously obtained by the numerical procedures 
using the kernel function (see, for example, ref. 20). 

With regard to boundary-value problems for finite wings in compress- 
ible flow, it appears that the procedure of separation of variables could 
be a feasible approach only for wings of very special plan forms such as 
a circle or an ellipse. In any case the development of the appropriate 
mathematical functions for a particular plan form would become highly 
involved. On the other hand it appears that approximate procedures 
similar to those used for two-dimensional wings might afford an approach 
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to solutions of these problems which though laborious might be handled 
by routine numerical methods. 

The kernel function of the integral equation relating pressure and 
downwash for the three-dimensional case appears as an improper integral 
that can be reduced and tabulated. Furthermore, its singularities can 
be isolated and determined so that use, similar to that made in the two- 
dimensional case, could be made of tabulated values of this function to 
obtain, by numerical procedures, aerodynamic forces for finite wings. 

The purpose of this paper is, therefore, to treat and discuss this 
kernel function. 


SYMBOIS 


h 0 ( 2) ,Hi( 2 ) 


io^i 


* 0*1 


K(x 0 ,y 0 ) 


K’(x 0 ,y 0 ) 

k 

IO,Ll 

L(l,n) 

2 

M 

P 


velocity of sound 

Hankel functions of second kind of zero and first order, 
respectively 

modified Bessel functions of first kind of zero and first 
order, respectively 

Bessel function of first kind of zero order 

modified Bessel functions of second kind of zero and first 
order, respectively 

kernel function of integral equation 
singular part of K(x 0 ,y 0 ) 

reduced-frequency parameter, Zcu/V 

modified Struve function of zero and first order, 
respectively 

unknown lift distribution 
reference length 
Mach number, V/c 
pressure 
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r 


- efa 2 + 



S 

t 

V 

w(x,y) 




region of xy-plane occupied by wing 
time 

forward velocity of wing 

amplitude function of prescribed downwash, 
v(x,y,t) = e iayt w(x,y) 

Cartesian coordinates 


x Q = x - | 


y 0 = y - n 

0 = {1 - M 2 

7 Euler's constant 



P 

t 

P 


03 



velocity potential 
acceleration potential 
fluid density 

circular frequency of oscillation 


ANALYSIS 

Integral Equation and Original Form of Kernel Function 


The main purpose of this analysis is to treat the kernel function of 
an integral equation that relates a known or prescribed downwash distri- 
bution to an unknown lift distribution for a harmonically oscillating 
finite wing in compressible subsonic flow. The integral equation referred 
to can be obtained by employing the Prandtl acceleration potential to 
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treat linearized boundary-value problems for oscillating finite wings by 
means of doublet distributions. Derivation of this integral equation 
from the linearized boundary -value problem for a wing is a preliminary 
task that has been done elsewhere (see, for example, ref. 21), but is 
reproduced herein as an appendix for the sake of completeness. 

In keeping with the concepts of linear theory, the wing is con- 
sidered a plane impenetrable surface S which lies nearly in the 
xy-plane as indicated in sketch l: 


z 



Sketch 1 


The x, y,z coordinate system and the surface S are assumed to move in 
the negative x-direction at a uniform velocity V. 

In terms of these coordinates the integral equation may be formally 
written as 

w(x,y) = i JJ 'l(M) K(x 0 ,y 0 ) d| drj 
S 


( 1 ) 
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where w(x,y) is the amplitude function of the prescribed downwash, 
K(x Q ,y 0 ) = K(x - l, y - ri) is the kernel function and physically repre- 
sents the contribution to downwash at a field point (x,y) due to a 
pulsating pressure doublet of unit strength located at any point 
and L(|,T|) is the unknown lift distribution or local doublet strength. 

The kernel function may be mathematically defined by the following 
improper integral expression (see eq. (A12), appendix A): 


K (Wo) = llin 

z — >0 


dz 2 


iasc 0 
e" V 



iao (a-m/a 2 + p 2 y 0 2 + P 2 z 2 ) 
\/a 2 + P^q 2 + P 2 z 2 


d A 


( 2 ) 


where M is Mach number, 3 = ]jl - M 2 , ciS = co/Vp 2 , a) is the circular 
frequency of oscillation, V is the velocity, and A is the variable of 
integration. Evaluation of this integral constitutes a main difficulty 
in obtaining aerodynamic coefficients for oscillating finite wings in 
compressible flow. The present analysis is therefore devoted to reducing 
it to a form that can be accurately evaluated by numerical procedures 
combined with the use of tables of certain tabulated functions. The form 
and order of all its singularities are determined and an expression for 
the kernel function is derived in which the singularities are isolated. 


Reduction of the Kernel Function 

In considering the reduction of the kernel function K(x 0 ,y 0 ), the 

integral involved can, for convenience, be written as the sum of two 
integrals, namely 

px n ianA-M A^+r 2 ) 
dA + / 2 dA 

J ° A 2 + r 2 

( 3 ) 

2 i^2 

V ( F 1 +F 2) w 

oz 


^x n ia) 
0 e 


.auA-M^+r 2 ) 


-icovA+Ml/A 2 +r 2 ) 


\/a 2 + r 2 


dA = 


0 Ja 2 + r 2 


Therefore, 


icux. 


K ( x o^o) = ljjn 


z — >0 


hi'* (F) = 


lim 
z — >0 
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where 


poo -ico (?v»-M \/a 2 + ) 

'W. ffT? " 


(5) 


and 


f 2 = 


p x o e io)(A-Ml/A 2 +r 2 ) 

0 W77 


dA 


(6) 


and where r = 3 + z^. 


The integrals F-j_ and F 2 are treated separately in succeeding 

sections. The final forms are given in equations (15) and (19), 
respectively. 

Evaluation of F^_. - The integral F^ can he converted to a form 
that can he more easily handled hy writing 


■i-r 


-ia)A e 


-i<joMVA 2 +r 2 


0 hfT? 


dA 


and introducing the following relation (see p. kl6 of ref. 22) 


-io5MlfA 2 +r 2 


\/a 2 + r 2 


= f J 0 (TA) 

^ 0 


-r^-M 2 ^ 2 

e 

I/t 2 - M 2 aj 2 


T dT 


'too 


_ J 0 (TA) ^ = ■ 

^T 2 - M 2 ” 2 


\/t^— M^ do 2 pMio -ir^M 2 aj 2 -T 2 

T dT - i / J n (TA) t dT 

J 0 fe 2 - T 2 


(JO 


( 7 ) 


In the first integral of these last two integrals make the subsitution 

^T 2 - M 2 ^ 2 = T 
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and in the second integral make the substitution 


'/m 2 !) 2 - T 2 = 


T 


Then 


e -iooM^A 2 +r 2 

!W77~ 





dT - 


i e”^ rT Jq(a^M 2 o } 2 - r 2 )dT 


(8) 


(It is of interest to note, in the expression on the left of eq. (8), 
that A and r appear in the same manner. The roles of these two quan- 
tities could therefore be interchanged in the expression on the right. ) 

With use of equation (8), the equation for F;j_ can be written as 


Fi = 




M 2 a> 2 






( 9 ) 


Changing the order of integration in each integral (which is a legitimate 
step because the integrands involved satisfy the continuity conditions 
required for such operations) leads to the following expression for F^ 
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The integrals within the brackets in equation (10) may be evaluated from 
tables of Fourier or Iaplace transforms as (see, for example, pair no. 55 
of appendix III of ref . 23 ) 


so that 



(id 


The first integral in equation (ll) can be written as 




-rT 


dT 


i 


T 2 - p 2 tD 2 



e 


~rr 




dT 


or 



-rT 


^/t 2 - p 2 <£ 2 


dT 


-Post cosh 9 



e -P<5r cos 0 


d9 


(11a) 


The first integral on the right of equation (lla) is given on page l8l 
of reference 22 as 


,-Pair cosh 9 


0 


d9 = Ko(pajr) 


where Kq is the modified Bessel function of the second kind of zero 

order. The second integral on the right of equation (lla) is given on 
page 338 of reference 22 as 

e-* 5 ^ 059 d8 --fjjotlfir) - I^teaarTl 
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where Iq is the modified Bessel function of the first kind of zero 
order and Iq is the modified Struve function of zero order. Then, 
the first integral of equation (ll) can be written as 



( 12 ) 


Note that the end result indicated in equation (12) is independent of 
Mach number. The second integral in equation (ll) may be written in 
another form as 



This integral has not been reduced to closed form; however, it is non- 
singular and can be readily handled by numerical methods. 

Combining equations (12) and (13) gives the following expression 
for F^: 



By performing the differentiations indicated in equation (4), there is 
obtained for the first part of equation (4) the following expression: 
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11m 




►0 dz £ 




iMp|yo| 

i e" P V 

3 



i/r 


2 

+ e 


— y 

v 4X0 



(15) 


All terms of this expression other than the integral may he evaluated 
at small intervals of y Q from existing tables, except at y Q = 0 where 

the function is singular. The integral is well-behaved and can be accu- 
rately evaluated by numerical or approximate procedures. The type and 
order of the singularities at y Q = 0 are discussed in a later section. 

Evaluation of F 2 . - In order to reduce the integral F 2 , equation (6), 
it is convenient to make the substitution 

A = r sinh 0 (l6) 


so that 


F o = 


\i 


-vsinh“-^!r 


'0 


e icur (sinh 0 -M cosh 0 )^g 


(IT) 


Noting that z appears only in r and performing the differentiations 
indicated in equation (4 ) yields 


= iso 

dz 2 /z=0 i y o| 


n . . -1 x 0 

Elnh w 


(sinh 0 - M cosh e)e^h°l < slnh 6 ‘ Mc0Sh O) d0 , £ ■— A* 

y 0 2/xo2+pSy 0 2 


-1_X0_ 

P|y Q | 


*° -is(x 0 -M^pgy 0 g) _Jm_ r e 2 coshe .(cosh0-MsirP,0)]e ii53 l yo l (sinhe - MCOShe) d9 

MP|y °' 


iS (x 0 pT/tT ) 







e iJi(x 0 -Ml(x 0 2 -p 2 y 0 2 )_ e -icBM0|y o | 


1 Phoi 

tap I 

cosh e e 1<; “ p l y ol (slnhe - Mcosh e 4e ► 


M M / 


- — 1 

1 J 

0 


(18) 
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or, by reverting completely to Cartesian coordinates through equation (16), 
there is obtained 




e iiu(x 0 -M|(xo 2 +P 2 yo 2 ) 
y 0 2 l/ x o 2 + P 2 y 0 2 


+ 




^ id) (x Q -M |/x q 1 2 + p ^y O 2 ) 


e -icuMp|y 0 



i£5 

e 




(19) 


This expression vanishes, as it should, for x 0 = 0 and like that 
in equation (15) has singularities at y Q = 0 which, also, will be 
handled in a later section. The integral that remains, like the integral 
remaining in equation (15), is nonsingular and simple in form and can be 
readily evaluated by numerical procedures. 


Expression for the kernel in terms of nondimens ional length variables . - 
Equations (15) and (19) can now be combined to give a reduced form of the 
kernel function K(x 0 ,y 0 ). However, in application, the variables x Q 
and y Q are employed, for convenience, in nondimens ional form. This is 
accomplished by considering these variables in a new sense to mean that 
they have been referred to some chosen length l and by introducing the 
reduced-frequency parameter k = Zcd/V. The variables will be used in this 
new sense throughout the remainder of the paper. The kernel can be written 
in terms of these nondimens ional variables as 


K(x 0 ,y 0 ) 


-ikx Q 

e 



F 2) 


z=0 






, , iMk ko 

1Mk |y 0 | + p — p 


m(*y 0 y 


r" /B ,C 2 -ik|y„|T , Mkx c * J(kx 0 )2 * S 2(ky 0 )2 £ K’" 

/ VI + T e dr - ■ ■ - e + 

° M (ky 0 ) 2 /(kx 0 ) 2 + P 2 (ky 0 ) 2 

1 p o 

M(ky o ) 2 ^0 
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Note that this expression for K(x 0 ,y 0 ) can be considered as a 
function of only three parameters, namely, k|y 0 |, kx Q , and M. To 
be more specific, the first two terms are functions only of k|y Q |; the 
next two terms are functions of k|y Q | and Mj and the last two terms 
are functions of k|y Q |, kx Q , and M. 

Equation (20 ) constitutes the principal result of this paper. Some 
partial checks as to its correctness are: (l) For k = 0, it reduces, 

as discussed subsequently, to the downwash of a pressure doublet in 
steady flow and (2) an integration with regard to the y-direction between 
the limits -» to +°o yields Possio^ result for the two-dimensional case. 
This integration is carried out in appendix B. Other special forms of 
the kernel function for M = 1, M = 0, and k = 0 are derived in subse- 
quent sections. In the section immediately following, the orders and 
types of the singularities of the kernel function are discussed. 


Discussion of the Singularities of the Kernel Function 


As previously indicated, the kernel function becomes singular or 
indeterminate at y Q = 0. The 


forms that the kernel function 
takes when it becomes singular 
are of particular importance in 
applications to lifting surface 
theory. It is therefore desir- 
able to extract and treat the 
s ingular i t ie s separately . 

This extraction can be 
conveniently made by con- 
sidering the value of K(x Q ,y 0 ), 
equation (20 ), at points on the 
semicircumference of a small 
ellipse (see sketch), the polar 


0 = -rt /2 


0 = 0 

— ► 


|0 = rt /2 

x Q Sketch 2 



equation of which may be written as 


x Q = e sin 0 

y Q = f cos 9 


(21) 


where, because of the symmetry of K(x 0 ,y 0 ) with respect to y Q , only 

the limits -n/2 $ 0 ^ rt/2 need he examined. Note that in these equa- 
tions values of 0 in the range -rt/2 ^ 0 < 0 correspond to field 
points ahead of or upstream from the doublet position and values of 0 
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in the range 0 < 0 £ «/2, to field points "behind or downstream from the 
doublet position. In particular, 9 = it/2 corresponds to points directly 
behind or in the wake of the doublet. 

After substituting these expressions for x© and y Q into equa- 
tion (20), the results may be written as 


K(e,0) = g 2 e lk€ 1 ke cos 9 K /ke cos 8 \ _ 

Z2 e 2cos20 I 3 \ 3 / 


iitke cos 9 

23 


/ke cos 9 \ T /ke cos 9 \ 

L Hr-j ' Ll (— — ) 


ikMe cos 9 

P 

M 


ike (sin 9-M) 

s 3 2 

M 


ikMe cos 9 


ike cos 9 p2 


sin 9 e 


ike (sin 9 -M) 

3 2 


. r'* £7? e-(^) T a T + 


P u 0 


pe sin 0 ^|(a-M^ 2 +€ 2 cos 2 9) 

/ e 3 2 dA 

mJ 0 


(22) 


With the use of the following series expressions for Kj_(z ) and 

[l x (z) - L^zjJ (which can be obtained from ref. 22 - for K^, see 
p. 80; for I 1 , see p. 77; and for L]_, see p. 329)'- 


%(*).(, + io g f)(f + ^ + fL + . ■ •) 


1 .(* + 54 + 25L + 

Z \k 6k 1152 


(23) 
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where 7 is Euler's constant (7 = 0.5772157), and 



z 2z 2 Z 3 2z^ z 5 
2 " 3n~ + 16 " 45* + 384 


(24) 


it is found that for vanishingly small values of e the limiting value 
of the expression for K(e,9) in equation (22) is for M < 1 


K(e,0) 


-ik€ sin 9 
e 





ke ( 1 - sin 0 ) 
2(1 - M) 



where 0(c n ) represents terms of order e 11 for n ^ 1. Expressed in 
terms of x 0 and y Q , equation (25) "becomes 





Examination of equation (25) shows that the kernel function K(e,9) has 
singularities with respect to e = ^/x Q 2 + as follows: 
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where, from equation ( 25 ), 


fl(0) = 


1 - sin 0 


(3 2 (1 + sin 0) 
cos^9 


> 


. . k(l - sin 0) 

f 2 (0) = log 7 -= log 


k cos 2 0 


2(1 - M) 2(1 - M)(l + sin 0) 


( 28 ) 


Although being of no particular significance in applications, it is 
of interest to note that the quantities f ^ and f 2 each have minimum 

values (lf-,1 =£ and |f 2 | , = log — at 0 = -n/2, which 

\' X| min 2 1 'min 1 - M/ 

corresponds to points directly ahead of the doublet positionj and, as 0 
increases from -jt/2 to +jt/2, the values of these quantities continu- 
ously increase from these minimum values to infinite quantities as 
follows : 


f l(l) 1 = lim 

5 — >0 

f3 2 1 + sing - 6) 

= lim 
6 — >0 

20 2 

cos 2 (— - 5^ 

8 2 


\2 ) 




> (29) 


»(! - u- 

6— >0 


k cos 


2 (i - 6 ) 


log 


2(1 - M) 


1 + sin 


(I - 5 )] 


= lim 
5 — ^0 


log 


kS‘ 


4(1 - M) 


Thus K(x 0 ,y 0 ) is singular for 0 = rt/2 even when the distance e from 

the doublet is not necessarily of zero order. This implies that the 
doublet produces a wake of discontinuous downwash that extends downstream 
from the doublet position to infinity. 

With knowledge of the singularities involved in the kernel func- 
tion K(x 0 ,y 0 ), an expression can be written in which the kernel is 
separated into a singular part and a nonsingular part (as was done by 
Schwarz, ref. 14, for the two-dimensional case) as follows 


K ( x o^o)= [K(x 0 ,y 0 ) - k’(wo 3 + K '( x o>y 0 ) 


(30) 
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where K(x 0 ,y 0 ) is defined in equation (20 ) or (22) and 



or in terms of e and 0, introduced by equations (21), 


-ike sin 9 


ik k 2 .. ke(l - sin 9) 
+ log — i - 

e 2 2(1 - M) 


The term [K(x 0 ,y o ) - K'(x 0 ,y 0 )J in equation (30) is a continuous func- 
tion for all values of k, x 0 , and y Q and for values of M in the 
range 0 <j M 1 . The term K'(x 0 ,y 0 ) is discontinuous at the doublet 
position (x 0 = 0, y 0 = 0) and at all points in the wake (x 0 > 0, 
y Q = 0 ). It is to be noted, however, that each term of K'(x 0 ,y 0 ) is 
integrable with respect to y Q or with respect to T] = y - y Q , a fact 
that may be useful in some numerical methods. 


Treatment of the Sonic Case 

Because of its special nature, the borderline case, M = 1, between 
subsonic and supersonic flow deserves and requires separate treatment. 

As M — > 1, the expression for the kernel function given in equa- 
tion (20) becomes indeterminate. It is possible, however, to obtain con- 
ditional limiting values for the kernel by considering the integral F, 
equation (4), and breaking it into two integrals, Fj_ and F2, as was 

done for the general case . 

With regard to F-j_, its limiting value and the value of its deriva- 
tives with respect to z at z = 0 can be shown to be zero as M — > 1. 
From the form of F;j_ given by equation (14), 
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But since (see ref. 22, p. 1J2 ) 

- f °° S = dr = -Kq(C ) (34 ) 

J 0 Vi+ t 2 

and (see ref. 22, p. 332) 



it may be concluded from equation (33) that 


lim F, = lim I = 1 = 0 

M — >1 M—>l\dz 2 / 


(36) 
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The total contribution to K(x 0 ,y 0 ) at M = 1, therefore, arises from 
the limit of F 2 , equation (6), as M — > 1. The limiting form of F 2 
may be written in terms of nondimens ional coordinates as 



In approaching the limit M = 1 (from the subsonic side) in equation (37), 
it is convenient to replace M by 


M = 1 - e 


where e is infinitesimally small so that 

p 2 = (1 - M)(l + M) = e (2 - € ) « 2e 
With this approximation, equation (37 ) may be written as 




(for x 0 > 0) (38) 


From physical considerations, the right side of equation (38) is to be 
considered zero for x 0 = 0. This is in keeping with results that would 
be obtained if the limit under consideration were sought from theory of 
supersonic flow, M > 1. 
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The integral in equation ( 38 ) cannot be completely expressed in 
terms of known functions. Furthermore, since it is singular at its 
lower limit, further treatment is required to reduce it to a form such 
that its derivatives with respect to z can be numerically evaluated. 
For this purpose the integral may be written as two integrals, namely 



= F 2 ' + f 2 " 


( 39 ) 


where 


and 



( 40 ) 




X 


dX 


( 41 ) 


The limits of integration in equation (4o) are so chosen that the inte- 
gral in this equation can be reduced to a known form, by making the 
substitution 



Thus, 
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Equation (42) may "be written in terms of the integrals involved in F 1 
(see eqs. (34) and (35))> namely. 



Differentiating this result twice with respect to z and then setting 
z = 0 gives 





(44) 


Differentiating equation (4l) twice with respect to z and setting 
z = 0 gives 


( 


^2 


It 







(45) 


After performing an integration by parts and collecting terms, equa- 
tion (45 ) may be written as 








e + 


k 2 y 0 2 




(46) 


Equations (44 ) and (46 ) are combined to give 



Then, in 


accordance with equation (4), there is obtained for K(x 0 ,y 0 ) M=1 : 
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For x Q > 0, 


K(xo,y 0 ) = •£- e " ikx °«J 


M=1 Z 2 


k |^c 


k 2 y 0 2 k 2 y 0 2 


( k M) 

iri 

T _ ( b I nr 

2 k |y 0 |L 

- L i^ K |y< 

k 2 S'o 2 '| 


^kx 0 

kx 0 j 

* + 1 

e 


k 2 y n 2 



J k |y 0 | 




dA^ ( 47 a) 


and> for x 5 0, 


K(x 0 ,y 0 ) M=1 = 0 


(47b) 


The integral appearing in equation (47a) is finite and proper and can 
be evaluated by numerical procedures. 


Treatment of the Steady and Incompressible Cases 


It is of interest to consider the form of the kernel function given 
in equation (20 ) for some particular values of M and k. In the fol- 
lowing sections a discussion is given for the steady case (k = 0) and 
the incompressible case (M = 0). The two-dimensional case is handled 
in appendix B. 


Reduction of the kernel for the case of steady flow . - In order to 
obtain the reduction of the kernel for the case of steady flow, consider 
the expanded form given by equation (26). As k — > 0, there results 
the following expression 


K(x °' y ° ) l t »0 



*o 


y Q 2 4o 2 + P 2 yo 2 / 


(48) 


which represents the downwash of a pressure doublet for steady flow. 
This result serves as a partial check as to the correctness of the 
expression for K(x 0 ,y 0 ) given by equation (20 ). 


By replacing y Q in equation (48) by y - and integrating 
from -1 to 1 with respect to tj, there is obtained 



i 

+ 

o 

L x 

'x 0 2 + p2(y - l) 2 

OJ 

x 0 (y - i) 


x„ + 


]fx Q 2 + p 2 (y + l) 2 

*o(y + 1 ) 


(49) 
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where the symbol 


/ 


indicates that the principal value or finite part 


of the improper integral must be taken. (See, for example, ref. 24 for 
a discussion of finite parts of such integrals. ) This result corre- 
sponds to the downwash produced by a simple horseshoe vortex two units 
wide. An equivalent expression for incompressible flow is given, for 
example, in reference 25, where in contrast to the present notation, 
x 0 has been chosen as positive forward. 


Reduction of, the kernel for M = 0 . - In order to effect the reduc- 
tion of the kernel for the incompressible case, the expressions for Fj_, 

equation (15), and F2, equation (l8), will be examined for the limit 
M — > 0: 

From equation (15) 

M^o & = t { K 4 k M)- f E.i( k M) - + *} (5°) 

z — >0 


and from equation (l8) 


OjFo 

lim 

M— > 0 dz 2 
z — >0 


ik 


inh 


-1 * 0 


|y 0 | 


sinh 9 e 


ik|y c 


I sinh 9 


d9 - 


x 0 ikx Q 

— e 

y 0 2 / x o 2 + y 0 2 


(51) 


Integrating by parts yields 


lim cfop. 


M— >0 


ik 


z — >0 S Z 2 


ik k 
^0 


■f ° i/y 0 2 + * 


2 e lkX 4X - 


ikXf 




y 0 


(52) 
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By setting x Q = 0 in equation (53 ), a form is obtained which can be 

shown to agree with results derived by Kussner for the case M = 0, 
x Q = 0 (ref. (26)). 

Some Remarks on Evaluation of the Kernel Function 

In regard to evaluating the kernel function for specific values 
of x Q , y Q , M, and k, an approximation for the function 

[K(x 0 ,y 0 ) - K'(x 0 ,y 0 ij 

can be obtained by making use of the series expansions for Kp (eq. (23)) 
and for (ip - Lp) (eq. (24)) and expanding all other terms of K(x 0 ,y 0 ) 

(eqi (20)) into a power series in terms of k. After performing these 
expansions and collecting terms with respect to powers of k, there is 
obtained 
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-ikx_ 


Z 2 3 2 



M - - - 1°.. - p 2 


i'-o 2 * e a y c 


fr-a- 


in 3^ 


ik5 

6p 2 


3 (l - 3M 2 )x q 2 + (2 - 3M 2 )p 2 y 0 2 
2 M^x 0 _ 


1923^ 


fo 2 + 3 2 y 0 2 
(12MP 2 - 20MV + 15P 6 - 12p 6 7 )y 0 2 - 32m5x 0 2 + 


4 ( 3 ^ + 6m 2 - l)x 0 3 + 12P 2 (m^ + 2M 2 - l) Xo y 0 2 

'/ x o 2 + e^o 2 




lk^ 

3603 6 


( 15 ^ + 10M 2 - l)x 0 2 |/x 0 2 + @ 2 y 0 2 - 4m3(5 + M 2 )x q 3 + 


3PV 0 4 


l/ x o 2 + P 2 y 0 2 


- 12M 5 p 2 Xo y 0 2 - 53 ^( 3 ^ - l)y 0 2 |/x 0 2 + p 2 y 0 2 


( 54 ) 


For values of the parameters that satisfy the following inequalities 

^o 


3 


< 1 


-§(x 0 - M||x 0 2 + 32y Q 2) < 1 


(55) 


equation ( 54 ) yields results that are correct to within 2 percent. For 
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values of the variables x Q , y Q , M, and k that do not satisfy the 

inequalities (55), values of the function [i(x 0 ,y 0 ) - K'(x 0 ,y 0 )] are 

perhaps most easily and economically obtained by numerically evaluating 
the two integrals in equation ( 20 ) and making use of available tables to 

evaluate other terms of |k(x 0 , y Q ) - K'(x 0 ,y 0 7|. 

It is pointed out that extensive tables of the Bessel functions 

and I-l may be found in reference 27 . A table of the Struve function L]_ 

with second and fourth differences for interpolation purposes may be found 
in reference 28 . 


CONCLUDING REMARKS 


The main purpose of this paper was to present the kernel function of 
the integral equation relating the downwash to the lift distribution in 
a form that can be computed. This purpose has been achieved by the pres- 
entation of the kernel in a form given in equation (20 ). This equation 
has been converted to a form more suitable for calculation by isolating 
the singularities as shown in equations (30) and ( 31 ). The special case 
of M = 1 is given in equations (47). The forms of the kernel function 
for other limiting cases, namely k = 0 and M = 0, are given in equa- 
tions (48) and ( 53 ), respectively. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va. , September 18, 1953. 
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APPENDIX A 

DERIVATION OF THE INTEGRAL EQUATION THAT RELATES THE DOWNWASH 
AND LIFT FOR A FINITE WING BASED ON REFERENCE 21 


In keeping with the concepts of linear theory the wing is considered 
as a nearly plane, impenetrable surface. Let this surface S lie nearly 
in the xy-plane, as indicated in sketch 1 of the body of the paper, and 
let it and the x,y,z coordinate system to which it is referred be assumed 
to move at a uniform speed V in the negative x-direction. At the same 
time let each point of the wing be assumed to undergo harmonic transla- 
tions of small amplitude Zm(x,y,t) at circular frequency co and let c 
represent velocity of sound in the medium. 

The problem for an oscillating wing consists in solving the wave 
equation subject to certain boundary conditions. The wave equation in 
rectangular coordinates is 


d 2 t djjr djjr 

ax 2 + c * 2 + * z 2 




0 


(Al) 


The independent variable in equation (Al) is regarded herein as 
an acceleration potential; as such it is directly proportional to a per- 
turbation pressure field and is related to a velocity potential p as 
follows : 


t = 


¥ + y ¥ 

dt dx 


(A2) 


In order to complete the boundary- value problem for the wing, it is 

, , dp 

desirable to calculate the downwash w(x,y,z,t) = — associated with i|r. 

oz 

Assuming this downwash to be. harmonic with regard to time implies that 
both potentials p and i|r are harmonic with regard to time and can be 
written, therefore, as 


jtf(x,y,z,t) 

i|r(x,y,z,t) 


e la)t 0(x,y,z) 
e iayt V(x,y,z ) 




(A3) 
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With these expressions for <f> and i|r, equation (A2) becomes independent 
of time and reduces to an ordinary equation with one dependent variable, 
namely 


dp 

\|r = icD0 + V — 
dx 


(AM 


This equation can be integrated with respect to x to give 


iosx 




icoA 

t(A,y,z)e v dA 


(A5) 


where the lower limit of integration is chosen, for later convenience, 
so as to satisfy the condition that vanish as x — > 

The boundary -value problem for the wing may now be expressed mathe- 
matically as follows: Under the assumption of harmonic motion the dif- 

ferential equation, equation (Al), becomes 



8 \|r 8 i|f 

8^ + a? 


C 2V 


v ! + 





0 


(a6) 


In order to insure tangential flow at the wing surface, the potential 
must satisfy the downwash condition 


w (x,y ) = (^j = ( v ^ + io^ZmCxjy) (A7 ) 

z— 0 

where w and Zm are amplitudes of velocity and displacements, respec- 
tively, and are assumed to be known from the motion of the wing. At 
z = 0, the pressure 

P = -p(M z=0 (A8) 

must be zero at all points (x,y) off the wing. At all points on the wing 
\|r is allowed to be discontinuous and the value of p at a given point 
is determined by the magnitude of the discontinuity in \| f at the point. 
In the neighborhood of the trailing edge, p must go to zero, corre- 
sponding to the Kutta condition. 
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One other condition, that fi vanish far ahead of the wing, is 
inherently satisfied by the relation between j) and \|r given in equa- 
tion (A5). 

The potential at point (x,y,z) due to a harmonically pulsating 

doublet located in the xy -plane at (|,tj,0) that satisfies equation (A 6) 
is 


\|r 

Y o 



icn t + 
e — 


c0 2 c(B 2 


R' 


(A9) 


where 

R' = \/(x - O 2 + P 2 (y - ti) 2 + p 2 z 2 

and. the factor A is a strength and dimensionality factor that makes 
possible different uses and interpretations of the potential ty 0 . If 

is considered as an acceleration potential and substituted into equa- 
tion (A5 ), there is obtained a corresponding velocity potential fi Q 
which may be written as 


P a = A A e V 
P ° dz 


, / +J . V MA R \ 

ico(x-|) , t 7; + ~ to To] 

- 77 -^ pA=x-| \ V cp2 C f3 2 / 


R 


dA (A10) 


where 


E - \fk 2 + 3 2 (y - I)) 2 + 3 2 z 2 


% 

The downwash associated with i(r 0 may be written as 

dz 


_ a £_ P x o e tt(*-Mh 2 +r 2 ) 

dz Zz? J + r 2 


dA 


(All) 
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where x Q = x - |, cd = co/vp 2 , and r = p>/(y - q) 2 + z 2 . With the use 

of this equation and the concept of solving linear boundary -value problems 
by means of superposition of elementary solutions to the governing dif- 
ferential equation, the boundary -value problem under discussion can be 
written as an integral equation, namely 


rr ic0X ° v2 r x Q i*U^) 

w(x,y) = lim A / / )e v d£ dq — - / dA (A12) 

z— *0 dz VA 2 + r 2 

where S represents the surface of the wing and L(£,q) represents an 
unknown lift distribution or doublet strength on S. Equation (A12) may 
be seen to correspond essentially to equations (l) and (2). 

If the distribution function L(|,q) in equation (A12) is deter- 
mined in accordance with the boundary conditions discussed in the pre- 
ceding paragraph, equation (A12) can be considered as a complete solu- 
tion to the boundary-value problem for an oscillating finite wing in 
compressible flow. It is also to be noted that equation (A12) can be 
considered to represent a solution to the so-called "indirect" problem, 
that is, that of finding the downwash distribution associated with a 
given lift distribution. 
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APPENDIX B 


REDUCTION OF THE KERNEL FUNCTION FOR THREE-DIMENSIONAL 
FLOW TO THAT FOR TWO-DIMENSIONAL FLOW 


The purpose of this appendix is to show that integration of the 
kernel function K(x 0 ,y 0 ) from -~ to +<» with respect to t\ = y - y 0 
leads to a known result for two-dimensional flow. The kernel is first 
modified to a form that, for the present case, is easier to handle. 

Then, after performing an integration by parts on the modified kernel, 
the form of the kernel for the two-dimensional case is given (eq. (Bl8)). 
In addition, the special cases of M = 1 (eq. (B23)) and M = 0 
(eq. (B30)) are also shown. 

The integration under consideration with respect to T} is equiva- 
lent to an integration with respect to y Q , namely 


l 



K(x 0 ,y-r]) dTj 



K(x 0 ,y 0 ) dy Q 


(Bl) 


It is remarked in advance that since z has been made zero in the 
expression for K(x 0 ,y 0 ), equation (20), it is necessary to employ the 
concept of "finite parts of infinite integrals" when integrating this 
function across the singularities at y Q = 0. Use of this concept gives 

the same results that could be obtained by the more arduous task of 
performing the integrations before setting z equal to zero. 

Modification of the kernel . - In order to effect the desired modifi- 
cation of the expression for K(x 0 ,y 0 ) given by equation (20), consider 
the first integral of the expression, namely 


This integral can be written as 

a* ** r e - T ( 5+lk i y ‘>i) d T ♦ r e- lk ' y ° |T a T 

5 — >0 Jo J M/p 


'M/p 


e- lk|y °l T 


dr 


(B2) 


0 



52 


NACA TN 3131 


but according to page 351 of reference 22 



e -T(6+ik|y 0 |) dT 


2(6 


-i-j— jy ^i(8 + Ik |y c | ) - Yi(» + Ik |y c ( j] 


(B4) 


where is the unmodified Struve function of first order and Yj is 

the Bessel function of the second kind of first order. In the limit as 
6 — > 0 these expressions have the following values: 

For the first expression in the bracket (see ref. 22, p. 329) 


lim E ± (b + ik|y 0 |) = H 1 (ik|y 0 |) = -L 1 (k|y Q |) (B5) 

6— >0 

and for the second expression (see ref. 22, pp. 77 and j8) 

lim Y X (8 + ik|y 0 |) = - 1 ^ ( 1 > (ik | y 0 1 ) + iJi(ik|jr 0 |) 

6 — >0 

“ T K l(k|y C |) - I^klYol) (B6) 


where denotes the Hankel function of the first kind of first 

order. With the use of equations (B3 ) to (B6), expression (B2) can be 
written as 



J: 


1+ T £ 


e- lk l y °l T dx 


k 2 f \|l+T 2 e " ik ^ y °l T dr + 


k 

ho| 



k l y o|) + f-[ii( k M)-k(k|y o 0 


(B7) 
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Substituting this result into equation (20) of the text gives the modi- 
fied. form of K(x 0 ,y 0 ) sought, namely 


K(x 0 ,y 0 ) = 


-ikx n 


iMk|y 0 | 

e 0 

ik 

e 3 + 1 

Z 2 

p|y Q | 

%o 2 


iMk|y 0 | 


MXq + \fx Q 2 + 3 2 y Q 2 ■p ( x o ~ M l/ x o 2 + 3 2 y o 2 ) 

^o 2 / x o 2 + 3V 


k 


2 



+ T c 


- ik |y 0 l T 

e 





) 


dA 


(B8) 

Integration of modified kernel . - Since the expression for K(x 0 ,y 0 ) 
is symmetrical with respect to y Q , that is, K(x 0 ,-y 0 ) = K(x 0 ,+y 0 ), 
the integration under consideration can be expressed as 


if k(x q , | y 0 |) dy Q = 
-00 



(B9) 


where, on the right, the absolute -value signs on y Q can be dropped. 
After performing an integration by parts by letting 


e 




-ikx Q 

e 


1 


dv = 2 


Z 


v = -2 


i y Q 


(BIO) 
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and 
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This expression vanishes at its upper limit y Q = » and is singular at 
its lower limit y 0 = 0. However, by not making z — > 0 in the deri- 
vation of K(x 0 ,y 0 ) until after this stage is reached, this singular 
value is canceled by other terms that have otherwise been dropped. Thus, 
the expression (B13 ) may be considered to be zero, which is the value of 
its finite part. The integration under consideration is then reduced to 
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The terms of this expression are treated separately in the next 
three equations: 

First (see ref. 22, p. l8o) 


lkM 

x o 


ik 


0 2 + P 2 y, 2 


... 


-i^|x 0 |cosh9 


dy 0 = ^ eP 2 M cosh 6+1^ 


d0 


second. 



(B15) 


\to -ihy 0 T 


2k 


dy 


0 °^M/p \[i + t 2 


dT = 2k c 


dT 


M/p \/l + t 2 J 0 


-iky 0 T 
e dy 0 


= -2ik 


CO 

r dT 

Jm/p t {i7T 2 


= -2ik log 


1 + P 
M 


(Bl6) 


and third (see ref, 22, p. l8o) 



illx-M^+P^o 2 ) 

^ „* - f 

ikA 

“tt 2 ax. 

r* - p2 ... dv. 

\j-K 2 + P^q 2 J ( 

0 Jk 2 + P^q 2 

CVJ 1 
CVJ 1 
II 

X ° 

e^ 2 dA 

0= -i|^lA|coshe 
/ e P d0 

p J, 

D J 

ikA 

0 

- * lk2 
p 

r° eP 2 H 
J 0 

o (2) 0f |x| ) dx (B17) 
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Substituting the results in equations (BI5) to (B17) into equation (Bl4) 
gives 


l r K(x 0 ,y 0 ) dy Q = e 1 

^ — 00 L P 


-ikx^ 


±kx r 


<e 


iM l x o| w (2)/kM|x 


H-, 

x„ 1 


O 1 )- 


H (2)/kMjXof 

0 \ p 2 , 


+ 3 log i±-6 + 

rt M 


ikA 

ik / ° e p2 H~( 2 


J 0 


*0 


Cf 1 )" 


(B18) 


This result is a form of the expression for the kernel function of 
Possio's integral equation relating pressure and downwash for a two- 
dimensional oscillating wing in subsonic compressible flow. It checks 
the results given, for example, in reference 29- 

Reduction of kernel for M = 1 .- The kernel function for M = 1 may 
be written as (see eq. (47a) ) 


K ( x °' y o )„ =1 



The second integral appearing in this equation can be shown to cancel 
several of the terms so that the kernel becomes 



K(xo,y 0 ), 


(B20) 
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so that the kernel for the sonic case in two-dimensional flow may he 
written as 

2 


-ikx. 


ikx r 


l / K(x 0 ,y 0 ) M=1 dy Q - 


2e 


ik y Q g 

2x, 




d y 0 - 


kx r 


iA 


i / e d dA 

J 0 '■ 


ik 2 y 0 2 

r 00 e" 2A 

— d V 

-«* y Q 


(B21) 


Integrating equation (B21) hy parts with respect to y Q , retaining only 
finite parts of the integrated results, and making use of the relation 


4 f°° -ia 2 - 2 

1 J e 


t2 dT = 2i f 

J 0 


^-ia 2 -!- 2 d T = ^ yields 



ikxp J 
2e 2 i 

( iky o 2 

00 

iky 0 2 \ 

e _ikX ° 

| J < 

1 " 2x 0 
e 

1 y 0 

ik 

“v 

-00 

C« *) 







iA 

r^o T 
i / e dA 

Jo 


ik 2 y c 2 


1 2A 
e 

y 0 


Finally, the kernel for the sonic case in two-dimensional flow may he 
written as 

/ ikx Q ,[hx^ I*?, 2 

f 00 , s . _ 4.|«I " ikx o 

z / K(x 0 ,yo)^^ dy Q 2 


ik^o 2 


ik^ f°° " 2 k 3 I 1 

_ e dy Q y 


ik 


JL. e 2 - ikl/iT / ' * e 2 dA I (B 23 ) 


•It may he noted that the integrals in this equation are readily express- 
ible in terms of Fresnel integrals 

ox 


and 


C(x) = / cos | t 2 dt 

y 0 


S(x) = I sin | t 2 dt 

Jq 2 
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Reduction of kernel for M = 0 .- For M = 0 it is convenient to 
modify the kernel function before integrating with respect to y 0 . For 
this purpose use is made of the relation (see eq. (BT ) ): 


K i( k |y 0 |) - 2TriEi( k l y ol) - L l( k l y ol j] = k2 


JT7T 2 e- ik l y °' T dr 


-±sf (B2k) 


y Q 2j 0 


and the relation 


^ r° fTTrf d x . ii r° 

y 0 2 Jo y Q 2 J-x r 


^y 0 2 + * 2 e ~ ikA ( B 25 ) 


With these relations the expression for K(x 0 ,y Q , equation (53), can 
be written as 


K ( x o^oW 0 = 


x o -ikx~ ik|/x 0 2 + y Q 2 ikx' 

y 0 2 / x o 2 + y G 2 y ° 2 


y o 2y -x c 


/y Q 2 + A 2 e" ikA dA 


(B26) 


But 



(B27) 
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therefore. 


K(x 0 ,y 0 )j^_o 


-ikx Q n c 0 


-ikA 


2 dA 

■ x ° (y Q 2 + a2) 5/2 


(B28) 


Integrating with respect to y Q gives 


‘L K( ^U^’-f e ' lta °/ 0 I. 


-ikA 


0 x °(y 0 2 + A2 ) 


3/2 


dA dy 0 


2 e- lkX ° f e- lkX d* 

l 


dy 0 


-X r 


0 (y 0 2 - x 2 f /2 


2 -ikx 0 P“ -ikA 

7 6 / 

1 J -x 0 A 2 


o i Vy / P 00 -ikA p 00 ikA 

-!•' °(I.T a -I ^ 4X| (B29) 


A 2 


Integrating each integral in equation (B29) and retaining only finite 
parts yields 


r“ . . 2 -itao I e 

X / ^-( x o^yo)M=n d y 0 i e 1* x 


ikx^ 


2 -ikx 0 / e' lkX ° 

= 'T e 


-ik* r“e lkA 1 

“LV“-“/ fc Tr“ 


d 0 


_ 2k slnjk* d>v _ 


ik f COS kA d7v + k f 3ln k A d>v 

J ^ J x. 


Wk _ 1 + _i_ e ’ ikx o, 

Z 2jtkx o 2 * 


jci(kx 0 ) + i|si(kx 0 ) + ft (550) 
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where -Ci(kx 0 ) and Si(kx 0 ) denote, respectively, the "cosine integral" 
and "sine integral' 1 functions defined as follows: 


Ci (x ) 




The results in the braces of equation (B30) check with results given 
for this case in reference l4. 
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